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We develop an analytical device model for graphene bilayer field-effect transistors (GBL-FETs) 
with the back and top gates. The model is based on the Boltzmann equation for the electron 
transport and the Poisson equation in the weak nonlocality approximation for the potential in the 
GBL-FET channel. The potential distributions in the GBL-FET channel are found analytically. 
The source-drain current in GBL-FETs and their transconductance are expressed in terms of the 
geometrical parameters and applied voltages by analytical formulas in the most important limiting 
cases. These formulas explicitly account for the short-gate effect and the effect of drain-induced 
barrier lowering. The parameters characterizing the strength of these effects are derived. It is 
shown that the GBL-FET transconductance exhibits a pronounced maximum as a function of the 
top-gate voltage swing. The interplay of the short-gate effect and the electron collisions results in a 
nonmonotonic dependence of the transconductance on the top-gate length. 
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I. INTRODUCTION 



Unique properties of graphene layers, graphene 
nanoribbon arrays and graphene bilayers [l|-|3| as well 
as graphene nanomeshs jj make them promising for dif- 
ferent nanoelectronic device applications. The gapless 
energy spectrum of graphene layers allows to use them 
in terahertz and midinfrared detectors and lasers 
9]. However, the gapless energy spectrum of GLs is 
an obstacle for creating transistor digital circuits based 
on graphene field-effect transistors (G-FETs) due to 
relatively strong interband tunneling in the FET off- 
state [ljj, EH- The reinstatement of the energy gap 
in graphene-based structures like graphene nanoribbons, 
graphene nanomeshs, and graphene bilayers appears to 
be unavoidable to fabricate FETs with a sufficiently large 
on/off ratio. Recently, the device dc and ac characteris- 
tics of graphene nanoribbon and graphene bilayer FETs 
(which are referred to as GNR-FETs and GBL-FETs, 
respectively) were assessed both numerically and analyt- 
ically [IU|<|. The device characteristics of GNR-FETs 
operating in near ballistic and drift-diffusion regimes can 
be calculated analogously with those of nanowire- and 
carbon nanotube-FETs (see, for instance [20l - [22| and 
references therein). The GBL-FET characteristics can, 
in principle, be found using the same approaches as those 
realized previously for more customary FETs with a two- 
dimensional electron system in the channel [23M251 H3? 
-l3l|. However, some important features of GBL-FETs, 
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in particular, the dependence of both the electron density 
and the energy gap in different sections of the GBL-FET 
channel on the gate and drain voltages should be consid- 
ered [33 - l34j . as well as the "short-gate" effect and the 
drain- induced barrier lowering [29j . 

In this paper, we use a substantially generalized ver- 
sion of the GBL-FET analytical device model [13, OH to 
calculate the characteristics of GBL-FETs (the threshold 
voltages, current-voltage characteristics, and transcon- 
ductance) in different regimes and analyze the possibility 
of a significant improvement of the ultimate performance 
of these FETs by shortening of the gate and decreasing 
of the gate layer thickness. The device model under con- 
sideration, which presents the GBL-FET characteristics 
in closed analytical form, allows a simple and clear eval- 
uation of the ultimate performance of GBL-FETs and 
their comparison. The following effects are considered: 
(a) Dependences of the electron density and energy gap 
in different sections of the channel on the applied voltages 
and the inversion of the gated section charge; (b) Degen- 
eracy of the electron system, particularly, in the source 
and drain sections of the channel; (c) The short-gate ef- 
fect and the effect of drain-induced barrier lowering; (d) 
Electron scattering in the channel. 

Our model is based on the Boltzmann kinetic equation 
for the electron system in the GBL-FET and the Poisson 
equation in the weak nonlocality approximation [35l [36j j . 
The use of the latter allows us to find the potential dis- 
tributions along the channel in the most interesting cases 
and obtain the GBL-FET characteristics analytically. 

The paper is organized as follows. In Sec. II, the GBL- 
FET device model under consideration is presented and 
the features of GBL-FET operation are discussed. In 



2 



^V d >0<^ 



V b >0 



-^V t <0 
Source | 



L, 




Top gate 
Drain 



Graphene bilayer 



Back gate (substrate) 



FIG. 1: Schematic view of the GBL-FET structure. 



Sec. Ill, the main equations of the model are cited. The 
general formulas for the source-drain current and the 
GBL-FET transconductance simplified for the limiting 
cases (far below the threshold, near threshold, and at low 
top-gate voltages corresponding to the on-state) are also 
presented. In this section, the source-drain current and 
the GBL-FET transconductance are expressed in terms 
of the Fermi energy in the source and drain contacts and 
the height of the potential barrier in the channel. To find 
the barrier height, the Poisson equation is solved for dif- 
ferent limiting cases in Sec. IV. The obtained potential 
distributions are used for the derivation of the explicit 
formulas for source-drain current and the transconduc- 
tance as functions of the applied voltages and geometri- 
cal parameters. Section V deals with a brief discussion of 
some effects (role of the device geometry, electron scat- 
tering, charge inversion in the channel, and interband 
tunneling), which influence the GBL-FET characteris- 
tics. In Sec. VI, we draw the main conclusions. Some 
reference data related to the voltage dependences of the 
Fermi energy and the energy gap in different sections of 
the GBL-FET channel are singled out to the Appendix. 



II. DEVICE MODEL AND FEATURES OF 
OPERATION 

We consider a GBL-FET with the structure shown in 
Fig. 1. It is assumed that the back gate, which is pos- 
itively biased by the pertinent voltage Vb > 0, provides 
the formation of the electron channel in the GBL be- 
tween the Ohmic source and drain contacts. A relatively 
short top gate serves to control the source-drain current 
by forming the potential barrier (its height A m depends 
on the top gate voltage V* and other voltages) for the 
electrons propagating between the contacts. 

We shall assume that the GBL-FETs under the condi- 
tions when the electron systems in the source and drain 
sections are degenerate, i.e., Ej? > fc^T. This implies 
that the back gate voltage is sufficiently high to induce 
necessary electron density in the source and drain sec- 
tions. 



In the GBL-FET the energ y g ap is electrically induced 
by the back gate voltage (32h34| (see also [Hj]). Thus in 
GBL-FETs, the back gate plays the dual role: it provides 
the formation of the electron channel and the energy gap. 
Since the electric field component directed perpendicular 
to the GBL plane in the channel section below the top 
gate (gated section) is determined by both Vb and Vt, the 
energy gap can be different in different sections of the 
GBL channel: E g>s (source section), E g (gated section), 
and E g _d (drain section) [13, [3 • At sufficiently strong 
top-gate voltage (V t < V t h < 0, where V t h is the thresh- 
old voltage), the gated section becomes depleted. Since 
the energy gaps in GBLs are in reality not particularly 
wide, at further moderate increase in |V t |, the gated sec- 
tion of the channel becomes filled with holes (inversion 
of the charge in the gated section) if V t < V, n < Vth, 
where V, n is the inversion voltage. As a result, the GBL- 
FETs with moderate energy gap are characterized by the 
threshold and inversion voltages: Vth < and V% n < 0. 
The explicit formulas for Vth and Vi n shall be given in the 
following. The cases V t = Vth an d Vt = Vj„ correspond 
to the alignment of the Fermi level in the source section 
of the channel with the conduction band bottom and the 
valence band top, respectively, in the gated section. 



III. MAIN EQUATIONS OF THE MODEL 

Due to relatively high energy of optical phonons in 
graphene, the electron scattering in the GBL-FET chan- 
nel is primarily due to disorder and acoustic phonons. 
Considering such quasielastic scattering, the quasiclas- 
sical Boltzmann kinetic equation governing the steady 
state electron distribution function / p = f p (x) in the 
gated section of the channel can be presented as 



dfp 

dx 



dtp df p 
dx dp x 



d 2 qw(q)(f p+(1 - /p)<5(£p^ 



(1) 



Here e = |e| is the electron charge, e p — p 2 /2r 
0.04mo is the electron effective mass in GBL (mo is the 
bare electron mass), p = (p x ,p y ) is the electron momen- 
tum in the GBL plane (z — 0), w(q) is the probability of 
the electron scattering on disorder and acoustic phonons 
with the variation of the electron momentum by quantity 
q = (q x ,Qy), v x = p x , and axis x is directed in this plane 
(from the source contact to the drain contact, i.e., in the 
direction of the current). For simplicity, we disregard the 
effect of "Mexican hat" (see, for instance, Ref. j34|) and 
a deviation of the real energy spectrum in the GBL from 
the parabolic one (the latter can be marked in the source 
and drain sections with relatively high Fermi energies). 
One of the potential advantages of GBL-FETs is the pos- 
sibility of ballistic transport even if the top-gate length 
L t is not small. In such GBL-FETs, one can neglect the 
right-hand side term in Eq. (1). 

As in Ref. [l(| EE HH HH , we use the following equa- 
tion for the electric potential tp — ip(x) = ip{x, z)\ z —q in 
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the GBL plane: 
(W b + W t ) d 2 <p ip-V b <p- Vt 



dx 2 



w b 



Wi 



47r e , 



(2) 

Here, E_ and £+ are the electron and hole sheet densities 
in the channel, respectively, k is the dielectric constant of 
the layers between the GBL and the gates, and Wb and 
Wt are the thicknesses of this layers. In the following, we 
put W b = W f = W (except Sec. V). Equation (2) is a con- 
sequence of the two-dimensional Poisson equation for the 
electric potential tp( x i z ) m the GBL-FET gated section 
(-L t /2 < x < L t /2 and -W b < z <W t , where L t is the 
length of the top gate) in the weak nonlocality approxi- 
mation [35] ]. This equation provides the potential distri- 
butions, which can be obtained from the two-dimensional 
Poisson equation by expansion in powers of the parame- 
ter 6 = {Wl + Wf)/lb{Wb + W t )/C 2 = W 2 /15C 2 , where 
£ is the characteristic scale of the lateral inhomogeneities 
(in the x-direction) assuming that S -C 1, i.e., C is not too 
small. The lowest approximation in such an expansion 
leads to the Shockley's gradual channel approximation, 
in which the first term in the left side of Eq. (2) is ne- 
glected [13, 38]. The factor 1/3 appeared due to features 
of the Green function of the Laplace operator in the case 
of the geometry under consideration. 

The boundary conditions for Eqs. (1) and (2) are pre- 
sented as 



= /.. 



p m >0,x=-L t /2 



fd, P , (3) 



P x<0,x=L t /2 



functions with the Fermi energies Ef, s and £F,di which 
are determined by the back gate and drain voltages, Vb 
and Vd [13, EH (see also the Appendix) : 



eV b 



e F ,d ~ e(V b - V d ) 



(5) 



where b = a B /8W, a B — kh 2 /me 2 is the Bohr radius, 
and H is the reduced Planck constant. In the following, 
we shall assume that b <C 1, so that ef, s — beVb and 
ep.d — be(Vb — Vd). In particular, if = 4 nm (GBL on 
Si02) and W = 10 nm, one obtains b ~ 0.05. Due to a 
smallness of 6, we shall disregard a distinction between 
V d * and Vd because Vd — V d * ~ bVd -C Vd (as shown in the 
Appendix) . Restricting ourselves by the consideration of 
GBL-FETs operation at not too high drain voltages, we 
also neglect the difference in the Fermi energies in the 
source and drain sections, i.e., put sp t d — £f,s = £f- 

The source-drain dc current density (current per unit 
length in the direction perpendicular its flow) can be cal- 
culated using the following formulae: 



J 



4e 



? -h 2 



{2nh) 2 
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oo JO 



dp y / dp x v x (f p - /_ p ). (6) 



In this case, Eq. (1) with boundary conditions (3) yield 



= 0, 



x=-L t /2 



x=L t /2 



= V d + {e F ,d -£ F ,s)/e = V d * , fp-f-p 



(4) 

where / SiP and fd, P are the electron distribution func- 
tions in the source and drain sections of the channel. 
The functions / SjP and fd, P are the Fermi distribution 



Q(pl/2m + etp) - @{jp 2 x /2m + ecp- eV d ) 
1 + exp[(p 2 /2m + eip - e F )/k B T] 



where T is the temperature, kg is the Boltzmann con- 
stant, and 6(e) is the unity step function. Using Eqs. (6) 
and (7), we obtain 



J 



ir 2 h 2 



dp v / d£ 



Am \ 1 + exp[(p 2 /2m + £ - e F )/k B T] 1 + exp[(p2/ 2 m + £ - e F + eV d )/k B T] 



ek B T 



dp v < In 



exp 



e F - P 2 y/2m - A,- 
k^f 



-In 



exp 



£f - P 2 y /2m - A m - eV d 
k~BT 



(8) 



Equation (8) can be presented in the following form: 



J = Jo / dz< In 



-In 



exp(£ m - z 2 ) + 1 



exp(<5 m - V d - z 2 ) + 1 



Here (see, for instance, Ref. |23l |) 

2V2^e(k B T) 3 / 2 



Jo 



ir 2 h 2 



(10) 



(9) 



is the characteristic current density, and S m = (ep — 
A m ) / 'k B T ', and Ud = eVd/k B T are the normalized volt- 
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FIG. 2: Normalized source-drain current J/ Jo versus normal- 
ized drain voltage Ud at different values of normalized gate- 
voltage swing 5 m . 



age swing and drain voltage, respectively. At m = 
4 x 1CT 29 g and T = 300 K, J ~ 2.443 A/cm. 

Figure 2 shows the dependences of the source-drain 
current J normalized by the value Jo as a function of the 
Ud calculated using Eq. (9) for different values of S m . 

The GBL-FET transconductance g is defined as 



9 = 



dJ 
Wt 



(11) 



Equations (9) and (11) yield 



g = J / dz 
Jo 



exp(z 2 - 5 m ) + 1 



exp(z 2 - d m + Ud) + 1 



d5 m 
dV t 



(12) 



The obtained formulas for the source-drain current and 
transconductance can be simplified in the following lim- 
iting cases: 

A. High top-gate voltages. At high top-gate voltages, 
which correspond to the sub-threshold voltage range, the 
barrier height exceeds the Fermi energy (A m 3> £ F ), so 
that S m ^> 1. In this case (the electron system in the 
gated section is nondegeneratc), using Eqs. (9) and (12), 
we obtain 



J 



9 = 



e F 



— Jo CX \ k B T 



/7F Jo / £f 

exp 



2 k B T 



k B T 



1— exp 



eV d \ 
k B T ) 



, (13) 



dA r , 



dV t 
(14) 



B. Near threshold top-gate voltages. In this case, A m > 
e F , i.e., \S m \ < 1, Eqs. (9) and (12) yield 



J ° cJ- aA 



k B T 



dV t 



(18) 



at high drain voltages eVd ~> k B T (Ud 1)- Here, ( — 
/ °°dein[exp(-^) + 1] ~ 0.678, Ci = / °° de/[exp(£ 2 ) + 
1] ~ 0.536, and (2 - J °° d£exp(£ 2 )/[exp(£ 2 ) + l] 2 . In 
the limit e F = A m , Eqs. (13) and (14) provide the values 
J and g close to those obtained from Eqs. (17) and (18), 
which are rigorous in such a limit. 

C. Low top-gate voltages. At low top-gate voltages, 
A m < sf, from Eqs. (9) 



J~ -I [(e F -A m f/ 2 e(e F - A r , 



- (s F - A rn - eV d ) 3/2 e(e F - A m - eV d )]. (19) 
g^I [(e F -A m )^ 2 e(e F -A m ) 

.(e F -A m -eV d f/ 2 @{e F -A m -eV d )][-^y 



Here. 



2V2me 

7T 2 ft 2 ' 



with Jo = I {k B Tf' 2 . 

Using Eqs. (19) and (20), one obtains 

J — lo eV d ^(e F - A m ) ^ h eV d y/sF~(l - ■ 



l T eV d 

9 - o J o -7= 

2 Je^ 



at eVd -C ef — A m , and 



dA m 
' dV t 



J^-Amf^hoeT^-'-t 



(20) 



(21) 



(22) 



(23) 



(24) 



■5 



g ~ I ^(e F - A m ) ^~ 



dA m 



(25) 



at eV d > £_p - A m . 

The dependences shown in Fig. 2 describe implic- 
itly the dependences of J calculated using the universal 
Eq. (9) on the back-gate, top-gate, and drain voltages as 
well as on the geometrical parameters. Equations (13) - 
(25) provide these dependences in most interesting lim- 
its. However, to obtain the explicit formulas for J as well 
as for g, one needs to determine the dependences of the 
barrier height A m on all voltages and geometrical param- 
eters. Since the electron densities in the gated section in 
the limiting cases under consideration arc different, the 
screening abilities of the electron system in this section 
and the potential distributions are also different. The 
latter leads to different A m vs V t relations. 



IV. POTENTIAL DISTRIBUTIONS, 
SOURCE-DRAIN CURRENT, AND 
TRANSCONDUCTANCE 

To obtain the explicit dependences of the source-drain 
current and the transconductance on the gate voltages 
Vb and V t as well as on the drain voltage Vd, one needs to 
find the relationship between the barrier height A m and 
these voltages. This necessitates the calculations of the 
potential distribution in the channel. The latter can be 
found from Eq. (2) in an analytical form in the following 
limiting cases. 



A. High top-gate voltages - sub-threshold voltage 
range (A m > e F ). 

When the barrier height A m exceeds the Fermi energy 
£_f, the electron density is low in the gated section and, 
hence, one can disregard the contribution of the electron 
charge in this section. In such a limit, we arrive at the 
following equation for the potential: 



ftp 
dx 2 



Fo 
A 2 ' 



(26) 



where A = ^/2/SW and F = -(V b + V t )/2. Solving 
Eq. (26) considering boundary conditions (4), for the case 
of high top-gate voltages we obtain 



ip = F 



cosh(x/Ao) 
cosh(L t /2A ) 



1 



Vd 



sinh[(2a; + L t )/2A ] 
sinh(L t /A ) 



(27) 

Limiting our consideration by the GBL-FETs with not 
too short top gate (L t >• W), Eq. (27) can be presented 
as 



; 2A^) C ° Sh (x 



Equation (28) yields 

V VaJ 2?7o 2 V VoJ 2 Va 

(29) 

where 770 = exp(L t /2Ao)/2. Simultaneously for the posi- 
tion of the barrier top one obtains 



Ao 
2 



2V d 



V b + V t 



(30) 



The terms in the right-hand side of Eq. (29) containing 
parameter % reflect the effect of the top-gate geometry 
(finiteness of its length). This effect is weakened with 
increasing top barrier length L t . The effect of drain- 
induced barrier lowering in the case under consideration 
is described by the last term in the right-hand side of 
Eq. (29). 

Equation (29) yields (dA m /dV t ) = -(e/2)(l - t^ 1 )- 
Invoking Eqs. (13) and (14), we obtain 



J = — Jo uxp 



e(V t 



2k B T 



Vth) f 1 _i_ 



1 — exp 



eV d 
k B T 



exp 



eV d 



2r l0 k B T 



(31) 



Vtt eJ 
9-—T ex P 



4 k B T 



e(V t - V th ) 



2k B T 



1 — exp — 



eV d 
knT 



exp 



1 

Vo 

eVd 
2r) k B T 



(32) 



Here, V th = -[1 + 26/(1 - v^Wb * -(1 + 2b)V b . The 
rightmost factors in the right-hand sides of Eqs. (31) 
and (32), associated with the effect of drain-induced 
barrier lowering, lead to an increase in g with increas- 
ing Vd not only at eVd ~ k B T but at eVd ^> k B T: 
g oc exp(eVd/2T] k B T). One can see that in the range 
of the top-gate voltages under consideration, the GBL- 
FET transconductance exponentially decreases with in- 
creasing \V t + V b \ and 



Je 
2k B T 



<5o 



\pH eJ 
~kB~f' 



(33) 



where at T = 300 K the characteristic value of the 
transconductance go ~ 4330 mS/mm. 



exp 



(28) 



B. Near threshold top-gate voltages (A m > ef)- 

-Y 

I&tf.JeVd < k B T <C Ef, taking into account that the 
electron distribution is characterized by the equilibrium 
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Fermi distribution function, the electron density in the 
gated section can be presented in the following form: 



2m , 



Considering this, we reduce Eq. (2) to 

d 2 (f <p _ F 



dx 2 A 2 A 2 ' 



(34) 



(35) 



Here. 



A 



a B W 
12(1 + 26) 



so that A/A ~ Vb < 1. The solution of Eq. (35) with 
boundary condition (4) is given by 



ip = F 



cosh(a;/A) 



cosh(L t /2A) 



V„ 



sinh[(2a; + J L t )/2A] 
sinh(L t /A) 



~ -F 



1 - 2(X1,( -2iJ C0Sh U 



Since A < Ao, one obtains rj 3> f]o, and the terms in 
Eq. (37) containing parameter r\ can be disregarded. This 
implies that the effects of top-gate geometry and drain- 
induced barrier lowering are much weaker (negligible) in 
the case of the top-gate voltages in question in compari- 
son with the case of high top-gate voltages. 

Substituting A m from Eq. (37) into Eq. (16), for low 
drain voltages we arrive at 



J e 2 V d ( 2 



At relatively high drain voltages (eV d > £f ^> ksT), 
the electron charge in the source portion of the gated 
section (x < x m , where x m is the coordinate of the bar- 
rier top) is primarily determined by the electrons injected 
from the source. The electron injection from the drain at 
high drain voltages is insignificant. Hence, the electron 
charge in the drain portion of the gated section can be 
disregarded. In this case, Eq. (2) can be presented as 



dx 2 

at —L t /2 < x < xq, and 



A 2 



F 
A 2 



(40) 



From Eq. (36) we obtain 



A m ~eF[ 1-- 
V 



eVd 
2 V 



(36) 



d 2 ^ 
dx 2 



A 2 



A 2 



(41) 



at x < x < L t /2. 

At the point x = xq corresponding to the condition 
e<^|x=x + £ f = 0, the solutions of Eqs. (40) and (41) 
should be matched: 



[e F -eb(V b +V t )} 1 



-ebVt 1 



V J 2r, 

(37) 

where r\ — exp(L t /2A)/2, and the position of the barrier 
top is 



\ eV d 






sf 






= V 




) 2 V ' 


X — Xq—O 


X — Xq-\-0 



2 F 2b V b ' 



(38) 



dip 
dx 



x—x^—0 



dtp 
dx 



(42) 



Solving Eqs. (40) and (41) with conditions (4) and (42), 
we obtain the following formulas for the potential ip at 
—Lt/2 < x < x and x < x < L t /2 as well as an 
equation for xq- 



<p = F 



cosh(a;/A) 
cosh(L t /2A) 



£j_ F cosh(x /A) 
e cosh(L t /2A) 



F 



cosh(x/A ) 



cosh(L t /2A ) 



V d 



sinh[(x + i t /2)/A)] 
sinh[(a;o + L t /2)/A)] ' 

sinh[(a; + L t /2)/A ] 
sinh( J L t /A ) 



(43) 



£f cosh(xo/A ) 
e °cosh(i t /2A ) 



F + V d 



smh[(x + L t /2)/A ] 
sinh( J L t /A ) 



sinh[(a; - L t /2)/A ] 
sinh[(a;o - L t /2)/A ] 



(44) 
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Here, as above, Vth — — (1 + 26) Vj,. In particular, 
Eqs. (48) and (49) at V t = V th , yield J th ~ J Co- For a 
GBL-FET with L t = 40 nm, W = 10 nm, at T = 300 K, 
one obtains J t h — 1.656 A/cm and gth — 2167 mS/mm. 



C. Low top-gate voltages (A m < ef). 

At relatively low top-gate voltages when A m < £p, 
the electron system is degenerate not only in the source 
and drain sections but in the gate section as well. In 
this top-gate voltage range, the spati al var iation of the 
potential is characterized by A ~ Wy/2b/Z. As a result, 
for A m one obtains an equation similar to Eq. (37). Since 
A < Ao <C L t , the parameter determining the effect of the 
top-gate geometry and the effect of drain-induced barrier 
lowering is r\ = exp(L t /2A)/2 ^> rj . As a consequence, 
one can neglect the effects in question in the top-gate 
voltage range under consideration. In this case, ef — 
A m ~ bV t is As a result, one can arrive at 

J ~ 2 -I Q e^[b{V t V th )} 3/2 - [b(V t - V th V d ]^ (50) 



when V d < b(V t - V th ), 



J^il e^[b{V t -V th )f'\ 



when V d > b(V t -V th ), 

OA, 



dV t 



-be. 



(51) 



(52) 



Considering Eq. (51), at low top-gate voltages we obtain 



bel eV d Vbe 3 / 2 I V d 



when eVd <C £f - A m ~ be{V b + V t ), and 
9 < bel a ^ ~ 6 3 / 2 e 3 / 2 / \/V b = 



(53) 



(54) 



when eV d ^> Ef — A m ~ be(Vb + V t ). As follows 
from Eqs. (33) and (34), the transconductance is pro- 
portional to a small parameter 6 3 / 2 . This is because 
the effect of the top-gate potential is weakened due to 
a strong screening by the degenerate electron system in 
the gated section. As a result, the transconductance at 
low top-gate voltages is smaller than that at the top- 
gate voltage corresponding to the threshold. Assuming 
that b = 0.05 and V b = 5 (e F ^ 0.25 eV), we ob- 
tain g on ~ 1467 mS/mm. Comparing Eqs. (49) and 
(54), we find g on /gth oc 6 3 / 2 V 'eV b /k B T ~ 0.158 and 
9 1 9th < goal gth ^ 0.68. 

Since the source-drain current at high top-gate volt- 
ages decreases exponentially when — V t increases, the 
transconductance decreases as well. 

Figure 3 shows the barrier (conduction band) profile 
A = -eip in the GBL-FETs calculated using Eqs. (27), 
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FIG. 4: Source-drain current J versus drain voltage Vd at 
different values of the top-gate voltage swing Vt — Vth- 



(36) , and (44) for different applied voltages and top-gate 
lengths. As demonstrated, the barrier heigh naturally 
decreases with increasing V t — Vth- At the threshold 
(Vt = Vth), the barrier heigh is equal to the Fermi en- 
ergy (at b = 0.05 and V b = 5 V, ep — 0.25 eV). One can 
see that shortening of the top-gate leads to a marked de- 
crease in the barrier height (the short-gate effect). The 
source-drain current MS cl function of the drain voltage 
for different top-gate voltage swings is demonstrated in 
Fig. 4. The dependences corresponding to V t — V t h < 0, 
Vt — Vth — 0, and V t — Vth > were calculated using for- 
mulas from subsections A, B, and C, respectively. Fig- 
ure 5 shows that the transconductance as a function of 
the top-gate voltage swing exhibits a pronounced maxi- 
mum at Vt ~ Vth- This maximum is attributed to the fol- 
lowing. At high top-gate voltages, the effect of screening 
is insignificant due to low electron density in the chan- 
nel. As a result, the height of the barrier top is rather 
sensitive to the top-gate voltage variations. The source- 
drain current in this case is exponentially small, so that 
the transconductance is small. In contrast, at low top- 
gate voltages, the screening by the electrons in the chan- 
nel is effective, leading to a much weaker control of the 
barrier height by the top voltage [pay attention to pa- 
rameter b <C 1 in Eqs. (50) - (54)]. Despite, a strong 
source-drain current provides a moderate values of the 
transconductance. However, in the near threshold volt- 
age range, both the sensitivity of the barrier height to the 
top-gate volatage and the source-drain current are fairly 
large. 



V. DISCUSSION 

A. Role of geometrical parameters 

In the main part of the paper, we assumed that the 
thicknesses of the gate layers are equal to each other: 
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FIG. 7: The same as in Fig. 6 but for Wt = 5 nm. 



770 > which strongly depends on the top- gate length L t . 
This effect gives rise to a decrease in the transconduc- 
tance if L t becomes smaller. However, in the range of 
large L t the electron collisions can play a substantial 
role. This leads to the transconductance roll-off with 
increasing L t . The effect of drain- induced barrier low- 
ering is also characterized by parameters 77 and 770. As 
follows from Sec. IV, at sub-threshold, the effect of drain- 
induced barrier lowering results in the appearance of the 
factor exp(eV d /2ri k B T) [see Eqs. (31) and (32)]. This 
factor can provide a marked increase in J and g with in- 
creasing drain voltage in GBL-FETs with relatively short 
top gates. For example, at W = 40 nm, L t = 40 nm, 
T = 300 K, and V d = 0.25 V, this factor is about 2.37. 



FIG. 6: Maximum transconductance g versus top-gate length 
L t for W t = 10 nm. 

W b = W t = W. If W b ^ W t , the formulas obtained 
above should be slightly modified. In particular, near 
the threshold (dA m /dV t ) = -(e/(l + W t /W b ))(l-rfo 1 ), 
so that the transconductance instead of Eq. (49) is given 
by 

Joe Ci ( 1 l\ 

9 k B T(l + W t /Wb)Co\ VoJ 

Joe Ci (5) 

k B T(l + W t /W b )Co' V ' 

As follows from the comparison of Eqs. (49) and (55), 
changing Wt, in particular, from W t = WbtoW t = Wb/2, 
leads to an increase in the transconductance at the top- 
gate voltages near the threshold of 50 percent (compare 
the g versus V t — Vth dependences in Fig. 5 for Wb — 
Wt — 10 nm and Wb — 10 nm and Wt — 5 nm. 

As demonstrated above, shortening of the top gate 
can result in deterioration of the GBL-FET character- 
istics. This effect is characterized by parameters 77 and 



B. Effect of electron scattering 

As shown above, the potential distributions in the 
main part of the gated section are fairly flat. This implies 
that to determine the effect of electron scattering asso- 
ciated with disorder and acoustic phonons one can use 
Eq. (2) with the collisional term following the approach 
applied in Ref. [17]]. Considering here the case when 
the elastic scattering mechanisms under consideration are 
strong, so that they lead to an effective isotropization of 
the electron distribution, one can find that the values of 
the source-drain current and the transconductance ob- 
tained in the previous section for the ballistic transport 
should be multiplied by a collisio n factor C. This factor 
is equals to C'co — y27r k B T /m/ L t v, where v = mw/2 
is the collision frequency (we put w{q) = w = const. It 
characterizes the fraction of the electrons injected into 
the gated section and those reflected back due to the col- 
lisions. To obtain the GBL-FET characteristics with the 
top-gate lengths in a wide range (to follow the transi- 
tion from the ballistic electron transport to the collision- 
dominated transport), we use for the collision factor the 
following interpolation formula: 
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1 + L t /L scat ' 

where L scat — \j2-KksT jm/v is the characteristic scat- 
tering length. Figures 6 and 7 show the dependences 
of the transconductance maximum (approximately at 
Vt = V t h) on the top- gate length calculated for GBL- 
FETs with different Wt- The scattering length is as- 
sumed to be oo (ballistic transport), 500 nm, and 75 nm. 
At T = 300 K, this corresponds to the collision frequen- 
cies v = 0, v ~ 1.14 x 10 12 s _1 (the electron mobili- 
ties n ~ 1.75 x 10 5 cm 2 /V s) and v ~ 7.6 x 10 12 s _1 
(/i ~ 2.63 x 10 4 cm 2 /V s) , respectively). One can see 
that in the case of essential electron collisions, g ver- 
sus L t dependences exhibit pronounced maxima. This is 
attributed to an interplay of two effects: the short-gate 
effect (weakening of the barrier controllability by the top- 
gate voltage when L t decreases) and the effect of colli- 
sions, which reinforces when Lt becomes larger, (which 
decreases the current). As follows from Figs. 6 and 7, the 
electron collisions can lead to a dramatic decrease in the 
transconductance. 



C. Charge inversion in the gated section. 

At sufficiently high top-gate voltages Vt < Vt n when 
A m > £p + Eg, the top of the valence band in the gated 
section of the channel can be markedly populated by 
holes (inversion of the gated section charge) , so that the 
term in the right-hand side of Eq. (2) becomes negative. 
The latter inequality corresponds to the following value 
of the inversion voltage [see Eqs.(A4) and (A5)]: 

V in = -V b (l + 2b+^c V th (l + ^) • (57) 

The hole charge provides an effective screening of the 
transverse electric field in the gated section. This leads 
to weakening of the sensitivity of the barrier height and 
the source-drain current on the top gate voltage V t [H[ 
and, hence, to a decrease in the transconductance. This 
pattern is valid at the dc voltages or when the char- 
acteristic time of their variation is long in comparison 
with the characteristic times of the thermogeneration of 
holes and the tunneling between the channel side regions 
and the gated section [l5|, Il7| . In the situation when 
the hole recharging of the gated section of the channel 
is a relatively slow process, the ac transconductance at 
the frequencies higher than the characteristic recharging 
frequency can substantially exceed the dc transconduc- 
tance. 



D. Interband tunneling 

At elevated top-gate voltages (V t < Vm), the inter- 
band tunneling of electrons from the conduction band 



in the source to the valence band in the gated section 
as well as from the valence band in the gated section 
to the conduction band in the drain can be essential. 
At high drain voltages the latter tunneling processes can 
be particularly pronounced. This can result in elevated 
source-drain current despite rather high barrier. To limit 
the tunneling, the back-gate voltage, which mainly de- 
termines the energy gap, should be high enough. A de- 
crease in the gate layer thicknesses also promotes the 
tunneling suppression. For example, if W = (5 — 10) nm, 
Vb = — Vt = 5 V, the energy gap in the gated section of 
the channel E g ~ 0.17 - 0.34 eV [see Eq. (A4)]. Despite 
some attempts to calculate the tunneling currents in G- 
FETs and GBL-FETs (see, for instance, Refs. [Til fl9j>. 
the problem for GBL-FETs remains open. This is be- 
cause the spatial nonuniformity of the energy gap in the 
channel and its nonlinear spatial dependence (particu- 
larly near the drain edge) associated with the features 
of the potential distribution under the applied voltages. 
Generalizing Eq. (A4), we obtain E g — ed[Vb— tp(x)]/2W, 
hence the energy gap varies from E g = E 9tS ~ edVb jW at 
x = -Li/2 to E g = E g<s ~ ed(V b - V d /2W at x = L t /2. 
It reaches a maximum E g = d(eVb + A m )/2W at x = x m . 
One can see that E g ^ can be markedly smaller than 
Eg, especially at not too small drain voltages. Due to 
this, the deliberation of the tunneling in GBL-FETs re- 
quires sufficiently rigorous device model (which could in- 
clude the formulas for the potential distribution obtained 
above) and numerical approach. 



VI. CONCLUSIONS 

We demonstrated that the developed device model of 
allows to derive the GBL-FET characteristics: the po- 
tential distributions along the channel and the depen- 
dences of the source-drain current and the transconduc- 
tance on the applied voltages and the geometrical pa- 
rameters as closed- form analytical expressions. The key 
element of the model, which provide an opportunity to 
solve the problem analytically, is the use of the Pois- 
son equation in the weak nonlocality approximation. In 
particular, the model accounts for the effect of screen- 
ing of the transverse electric field by the electron charge 
in the channel, the short-gate effect, and the effect of 
drain-induced barrier lowering. The parameters 770 and 
77 characterizing the strength of these effects in the cases 
of essentially depleted channel and strong screening were 
expressed via the geometrical parameters and the Bohr 
radius. As shown, the GBL-FET transconductance ex- 
hibits a pronounced maximum as a function of the top- 
gate voltage swing. The interplay of the short-gate ef- 
fect and the electron collisions results in a nonmono- 
tonic dependence of the transconductance on the top- 
gate length. The obtained analytical formulas for the po- 
tential barrier height, the source-drain current, and the 
transconductance can be used for GBL-FET optimiza- 
tion by proper choice of the thicknesses of gate layers, 
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the top-gate length, and the bias voltages. 
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Appendix 

Disregarding the effect of "Mexican-hat" and the non- 
parabolicity of the electron energy spectrum, the density 
of states can be considered independent of the energy (in 
the energy range under consideration). Taking this into 
account, the electron Fermi energies and the energy gaps 
in the source and drain sections of the GBL-FET channel 
are, respectively, given by [l?], Ell (Wi, = W t = W) 



k B T 
(1 + &) 



In 



exp 



beV b 
k R T 



- 1 



£ F,c 



k B T 
(1 + b) 



■In 



exp 



( be(V b -V d ) 
V k B T 



E„ 



edV b 
~2W~' 



E„ 



ed(V b - V d ) 
2W 



(Al) 



(A2) 



(A3) 



Here a B = kh 2 /me 2 , b = a B /8W, and d < do, where 
do — 0.34 nm [30j is the spacing between the graphene 



layers in the GBL, while d stands for the effective spac- 
ing accounting for the screening of the transverse electric 
field by GBL (polarization effect). In portion of the gated 
section essentially occupied by electrons and its depleted 
portion, one obtains 



E g = 



ed(V h - V t ) 
2W 



ed (V b - V t ) 
2W 



(A4) 



respectively. Due to > c!, from Eqs. (Al) - (A4) one 
obtains £f <s > Ep,d > E g<s > E g>d . At V t < 0, E g can 
significantly exceed E 9:S and E g _ d . In the case of strong 
degeneracy of the electron system, Eqs. (Al) and (A2) 
yield 



e F , s — eV b 



a B /8W 



£F,d ^ e(V b - Vd) 



(1 + b) 



(A5) 



The quantity Efa is given by the same equations in 
which, however, V b is substituted by V b — V d . As a result, 



£F,s — £F,d 



beV d 
1 + b 



zV d = eV d + e F ,d - £f, s 



beV d , 



eV d 



1 



Since parameter b in reality is small, so that Ef,i 
beV d < eV d and V^ ~ V d we put e F , s = £fa - 
substitute VT by V d . 



(A6) 



(A7) 



- e F ,d ^ 
Ef and 
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